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ABSTRACT 

We have discovered strong gravitational lensing by the galaxy ESO325-G004, in images obtained 
with the Advanced Camera for Surveys on the Hubble Space Telescope. The lens galaxy is a boxy 
group-dominant elliptical at z = 0.0345, making this the closest known galaxy-scale lensing system. 
The lensed object is very blue (B — I c w 1.1), and forms two prominent arcs and a less extended third 
image. The Einstein radius is i?Ein = 1.9 kpc (~ 3arcsec on the sky, cf. 12arcsec effective radius of 
the lens galaxy). Assuming a high redshift for the source, the mass within i?Ein is 1.4 x 10 11 M@, 
and the mass-to-light ratio is 1.8(M/L)qj. The equivalent velocity dispersion is a\ ens = 310km s -1 , 
in excellent agreement with the measured stellar dispersion a v = 320 km s _1 . Modeling the lensing 
potential with a singular isothermal ellipse (SIE), we find close agreement with the light distribution. 
The best fit SIE model reproduces the ellipticity of the lens galaxy to ~10%, and its position angle 
within 1°. The model predicts the broad features of the arc geometry as observed; the unlensed 
magnitude of the source is estimated at I c ~ 23.75. We suggest that one in ^200 similarly-massive 
galaxies within z < 0.1 will exhibit such a luminous multiply-imaged source. 
Subject headings: galaxies: elliptical and lenticular, cD — gravitational lensing 
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1. INTRODUCTION 

Gravitational lensing of background sources can yield 
valuable information on the mass profiles of galaxies, 
groups and clusters. In contrast to other methods, lens- 
ing constraints are independent of assumptions about the 
dynamical or hydrodynamic state of tracer material (e.g. 
stars, galaxies, X-ray gas). On galaxy scales, lensing 
helps to lift the degeneracy between the potential and 
the orbital anisotropy which plagues dynamical mass es- 
timates. For distant galaxies, lensing constrains the mass 
enclosed at large radii, beyond the reach of stellar dy- 
namical studies (e.g. Treu & Koopmans, 2004). 

The optimal configuration for lensing is with the de- 
flecting potential at half the distance to the source. 
Hence lensing is usually observed for systems at inter- 
mediate redshift, z ~ 0.3, where the number of po- 
tential background sources is very large. At low red- 
shifts (z < 0.1), strong lensing systems are rare. On 
cluster scales, a number of lensed arcs have been dis- 
cussed (e.g. Blakeslee et al. 2001). The nearest known 
galaxy-scale strong-lensing system, prior to this Letter, 
was Q2237+0305 (Huchra et al. 1985), a four-image 
QSO lensed by a z = 0.039 spiral. Among galaxy lenses 
with extended arcs, the lowest known lens redshift is 
z = 0.205, for SDSS1402+6321 (Bolton et al. 2005). 

1 Based on observations made with the NASA/ESA Hubble 
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which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program 10429. 
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This system was discovered on the basis of anomalous 
emission lines in the lens galaxy spectrum. A more direct 
approach to finding extended arcs is, of course, through 
imaging observations. Blakeslee et al. (2004) have dis- 
cussed the excellent prospects for serendipitous discov- 
ery of strong lens galaxies with the Advanced Camera 
for Surveys (ACS) on the Hubble Space Telescope. 

In this Letter, we report ACS discovery of a new 
galaxy-scale lens, with multiple extended images, at 
z = 0.0345. To our knowledge, the elliptical ES0325- 
G004 is the nearest known strong lensing galaxy, and 
provides for the first time a low-redshift analogue of dis- 
tant galaxy-scale arc systems. 

We adopt the WMAP cosmological parameters, i.e. 
H = 71km s^Mpc" 1 , fi m =0.27, fi A =0.73 (Bennett et 
al. 2003). 

2. THE LENS GALAXY ESO325-G004 

The massive boxy elliptical galaxy ESO325-G004 
(13' , 43 m 33!20, -38°10 m 33"6) lies at the center of the 
poor cluster Abell S0740 (Abell, Corwin & Olowin 1989), 
and is probably the dominant galaxy of that group. The 
galaxy has radial velocity cz =10420 km s _1 in the Cos- 
mic Microwave Background frame (Smith et al. 2000), 
corresponding to an angular scale of 0.678 kpc arcsec -1 . 
The galactic extinction is E(B — V) = 0.06 (Schlegel, 
Finkbeiner & Davis 1998). The stellar velocity dispersion 
of ESO325-G004 was measured by Smith et al. (2000), 
within an aperture 3.8x3.0 arcsec 2 ; the average of their 
two measurements is cr v = 320 ± 7 km s~P Smith et al. 
(2001) report the effective radius as R e s = 12.5 arcsec 
(8.5 kpc) from R-band imaging. 

Understanding the environment of the lens galaxy 
can be critical for correct interpretation of lensing con- 
straints. Figure Q shows the projected galaxy density 
around ESO325-G004, and the distribution of published 
redshifts. Although not complete, the redshift data sug- 
gest that S0740 is distinct from the neighboring cluster 
Abell 3570 (^40arcmin away, and with mean redshift 
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Fig. 1. — Environment of the lensing galaxy. The main fig- 
ure shows contours of projected galaxy density from the 2MASS 
extended source catalog. Small points mark galaxies with mea- 
sured redshifts in the NASA Extragalactic Database, within 1.5° 
of ESO325-G004. Galaxies with redshifts within 500 km s" 1 of 
ESO325-G004 are shaded. The histogram above shows the same 
redshifts, with the same redshift interval highlighted. 



-1000km s" 1 larger). 

ESO325-G004 was observed with ACS in January 
2005, as part of our program measuring surface- 
brightness fluctuations in the Shapley Foreground re- 
gion. The total integration times were 18882 seconds 
in F814W (22 individual exposures) and 1101 seconds 
(three exposures) in F475W. The frames were combined 
using MULTIDRIZZLE in STSDAS to yield the stacks used 
in this study (Figure |5J). Inspection of the images re- 
vealed two long, narrow arcs at —3 arcsec separation from 
the galaxy center (Figure EJ. A third object, less obvi- 
ously distorted, is present at similar radius. The arcs are 
prominent even in the much shallower F475W exposure 
(Figure Et), due to favorable color contrast between the 
blue arcs and the red light of the foreground elliptical. 

For a preliminary photometric analysis of the lens 
galaxy, we applied iraf tasks based on the ellipse-fitting 
algorithm of Jedrzejewski (1987). The surface photome- 
try model includes the C4 Fourier coefficient, which is nec- 
essary to describe the strong boxiness of the isophotes. 
All other high-order Fourier terms were forced to zero. 
The luminosity profile follows a i? 1 / 4 law out to at 
least —larcmin (~ 5i? ff). Transformed to the John- 
son/Cousins system (following Sirianni et al. 2005), and 




Fig. 2.— The deep (18882 s) F814W image of ESO325-G004 
and the surrounding field. The white square indicates the central 
8x8 arcsec 2 region covered by Figure [3] 



correcting for extinction and fc-dimming, the color is 
B — I c w 2.5, typical for an old, metal rich stellar popu- 
lation. At the radius of the arcs, the isophotal position 
angle is 68°, and the ellipticity e = 1 — ^ w 0.25. 

Although the arcs are seen clearly in the original im- 
ages (Figure^), their visibility is enhanced by subtract- 
ing a model for the lens galaxy light. We have explored a 
number of approaches to this subtraction. The simplest 
method uses the ellipse+C4 fits described above, sub- 
tracting independent models for the F814W and F475W 
images. The resulting F814W residual image is shown in 
Figure Ob- This method leaves a strong spoke pattern; 
there is also a risk of introducing tangential residuals un- 
related to lensing. A more robust approach is to use the 
color contrast between the arcs and the galaxy, subtract- 
ing a scaled version of the F814W image from the F475W 
image. The color residual image (Figure|3]:) is limited by 
noise from the shallow F475W exposure, but is system- 
atically very clean; in particular the spoke pattern and 
many of the point sources are removed. 

In Section |3J we will show that the observed geometry 
of arcs 'A', 'B' and 'C is indeed consistent with a single 
background source lensed by an elliptical potential. We 
have determined approximate magnitudes and colors for 
the arcs, based on the residual images, as summarized in 
TableQ] In particular, note that the arcs have consistent 
colors, with B — I c = 1.10 ± 0.05 (extinction corrected). 
Such a blue color can be reproduced with star-forming 
spectral templates at z < 0.45, with [O II A3727A] con- 
tributing to the F475W flux. The arcs have irregular 
surface brightness distributions, with numerous 'knots', 
suggestive of star-forming regions. In 'B' we can identify 
at least seven knots, while 'A' shows a symmetric struc- 
ture as expected for a two-image arc, and 'C appears 
double. 
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TABLE 1 

Properties of the arcs 




ID 


radius a 


lcngth a 


FSUW h 


F475VK b 


(-B -/"c)corr 


A 


2.8 


3.4 


22.06 


22.93 


1.07 


13 


3.4 


1.5 


22.92 


23.88 


1.21 


C 


2.6 


0.9 


22.61 


23.45 


1.03 




a Radius and length (in arcseconds) of an annular segment bound- 
ing the /^F8i4W = 23 mag arcsec -1 isophote, as used for the pho- 
tometry. 

integrated magnitude and color, on the Vega system, within the 
annular segment. 
c Johnson/Cousins, corrected for galactic extinction. 
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Fig. 3. — Observations and models of lensing in ESO325-G004. 
Panel (a) shows the 1101 s F475W image, with color map optimized 
to show the arcs. In Panel (b), the 18882 s F814W image is shown 
after subtracting a smooth boxy-elliptical profile model. In (c), we 
show the color-subtracted image, with results from a lensing model 
for a singular isothermal ellipse (see text). 



3. MASS-TO-LIGHT RATIO FROM LENSING 

We have used simple mass models to test the assump- 
tion that all three images are generated by gravitational 
lensing, and to derive initial estimates for the mass en- 
closed within the arcs. To test parametrized forms for the 
lensing potential, we applied the 'curve-fitting' method 
as implemented within the GRAVLENS/lensmodel soft- 
ware (Keeton 2001). This method takes as input a set 
of points on each observed arc, and optimizes a lensing 
potential which maps points from each curve to counter- 
images on the other curves. The input curves were ob- 
tained by manually tracing arcs 'A', 'B' and 'C. To de- 
termine a well-defined Einstein radius for an elliptical 
potential, Rusin, Kochanek & Keeton (2003), have sug- 
gested a scheme based on fitting a singular isothermal 
sphere (SIS) plus external shear. The model has three 
parameters, the SIS Einstein radius, and the magnitude 
and position angle of the shear. This parameter space 
was explored using a combination of grid-search and sim- 
plex optimization methods. 

The best SIS+shear fit yields i?Ein ~ 2.85arcsec 
(1.93 kpc or 0.23i? o ff). From this we can determine 
the mass interior to i?Ein to be M^in ~ l-4(£) s /£)i s ) x 
10 11 M Q . Here D s is the angular-diameter distance to 
the source, and D\ s the angular-diameter distance be- 
tween lens and source. The enclosed mass measurement 
thus depends on the source redshift, which is as yet un- 
known (but see below). For a given lens, the source is 
most likely to be at high redshift, where the surface den- 
sity of potential sources is very large. In this section, 
we consider the limiting case of a distant source, such 
that D s /D\ s w 1. The observed magnitudes interior to 
i? E in are F8UW = 13.5 and FA75W = 15.6. Trans- 
forming these to Johnson/Cousins bandpasses following 
Sirianni et al. (2005), we have I c = 12.7 and B — 15.1 
(corrected for extinction and /c-dimming). Converting 
to luminosities we have LEin,i c = 8.0 x 10 10 Lqj c and 
-^Ein.B = 3.0 x 10 10 L Q: b, and the mass-to-light ratios 
at the Einstein radius are TEi n .i c = 1-8 (M/L)@j c and 
TEin,B = 4.7(M/L)0 j s. Such values are typical for stel- 
lar mass-to-light ratios in old populations. Thus the mass 
budget within the Einstein radius appears to be domi- 
nated by the observed stellar mass 6 . The SIS mass model 

6 Given independent information on the stellar population (e.g. 
from spectroscopic line indices), this result could yield an upper 
limit on the mass contribution from dark matter associated with 
ESO325-G004 and/or the surrounding cluster S0740. 



4 



Russell J. Smith et al. 



yields ctsis = 310km s -1 , in excellent agreement with the 
measured stellar velocity dispersion. 

Considering the morphology of the galaxy, and the ex- 
pectation that the mass within i?Ein will be dominated 
by stars, a more realistic model for the lensing potential 
is the Singular Isothermal Ellipse (SIE; e.g. Korman, 
Schneider & Bartelmann 1994). In fitting the SIE model, 
we force the center of the potential to align with the ob- 
served galaxy, but allow the mass normalization, position 
angle and ellipticity to vary. Formally, the best fitting 
model has ellipticity e = 0.28 with position angle 68°, 
very similar to the equivalent parameters for the galaxy 
light. Figure 0t shows the results of this best-fitting SIE 
model. The reconstructed source position is ~0.3 arcsec 
north of the lens center (small circle) , and likely strad- 
dles the inner 'astroid' caustic of the lens model (inner 
cusped curve). The figure shows the predicted image lo- 
cations for such a source. Qualitatively, the SIE potential 
reproduces the broad features of the observed geometry, 
and confirms image 'C as a counter-image of 'A' and 
'B'. The long arc 'A' comprises two images crossing the 
critical curve (outer ellipse). While broadly successful, 
it is also clear that the model does not match the data 
perfectly at a more detailed level. In particular, the SIE 
does not match the relative lengths of the three observed 
arcs. For 'C, moreover, the model predicts a position 
displaced towards 'B', relative to the observed location. 
To improve the reconstruction, more sophisticated mod- 
els could take into account the intrinsic morphology of 
the source, the boxiness of the lens galaxy, and a possible 
external shear term. 

4. DISCUSSION 

It is interesting to consider whether the lensing configu- 
ration observed in ESO325-G004 is intrinsically unusual, 
or whether many other nearby galaxies might reveal such 
bright arcs when observed in sufficient detail. The cross- 
section for four-image lensing under the SIE mass-model 
is ~1 arcsec 2 , i.e. roughly the area of the 'astroid' caus- 
tic on the source plane. The estimated magnification 
factors for the separated arcs 'B' and 'C are ~4, suggest- 
ing the unlensed magnitude of the source is I c ~ 23.75. 
To this limit, the integrated I-band galaxy counts are 
~ 10 5 deg -2 (Postman et al. 1998), so the probability 
that a given galaxy has such a bright background source 
aligned for quadruple imaging is ~0.6%. To estimate 
the total number of massive galaxies available to serve 
as low-redshift lenses, we use the 2MASS J-band lumi- 
nosity function of Cole et al. (2001). ESO325-G004 has 
Mj — 51og/i = —24.16, and thus a luminosity ~ 7Lj. 
Integrating the luminosity function, parametrized as a 
Schechter function, the space density of galaxies above 
7L*j is 2.8 x 10~ 5 Mpc~ 3 (for h = 0.71). This density 
implies ~ 10 4 galaxies, at least as luminous as ES0325- 
G004, within z = 0.1. Combining these estimates, we 
find the total expected number of 'similar' low-redshift 
lens systems to be ~60, with only ~3 at the distance of 
ESO325-G004 or closer. 

As noted above, the lens model normalization depends 
on the source redshift, through the factor D s /D\ s . We 
have re-examined the raw spectra obtained by Smith 
et al. (2000) at the Anglo-Australian Telescope, which 
cover the range 4925 — 5740 A, with slit intercepting Arc 
'C. No anomalous emission lines are seen at the ex- 



pected position of the arc, which at face value suggests 
z SIC > 0.54 and z src < 0.32 (absence of [O n A3727A]), 
but also z src > 0.15 (absence of [O in A5007A] and H/3). 
Over the redshift interval 0.15 — 0.32, the lensing mass 
correction factor D s /D\ s ranges from 1.30 to 1.13. It is 
however quite possible that the exposures were too short 
(600 sec), and the slit too wide (3 arcsec), to detect emis- 
sion from the arc against the high background of the lens 
galaxy. 

Finally, we note that our images show four extremely 
faint tangential arc candidates at greater separation 
from ESO325-G004 (radius 9.4 arcsec, at position an- 
gles -95°, 90°, -35°, 150°). The surface brightness of 
these features is very low (> 24.5 mag arcsec^ 2 in I). If 
confirmed, the outer arcs could provide additional con- 
straints on the mass distribution at larger radius, and 
on the relative contribution of the group potential to the 
lensing mass. 

5. CONCLUSIONS 

We have discovered strong gravitational lensing by the 
elliptical galaxy ESO325-G004, which at z = 0.0345 
is the nearest known galaxy-scale lens. The multiply- 
imaged background source appears to be a star-forming 
galaxy, with prominent substructure. If the source is 
very distant, the arc radius is consistent with a stellar- 
dominated mass-to-light ratio within radius ~ \R c a- An 
elliptical isothermal mass model recovers the position an- 
gle and ellipticity of the galaxy, independent of the ob- 
served luminosity. The best-fit mass scale is consistent 
with the measured stellar velocity dispersion. Future 
modeling of the system should incorporate structural in- 
formation from the clumpy arc morphology. We intend 
to obtain integral-field spectroscopy for ESO325-G004, 
with 8m-class telescopes, to secure the source redshift 
and enclosed mass estimate. Additionally, these obser- 
vations will yield extended stellar dynamics for the lens 
galaxy, and enable high-contrast imaging of the arcs by 
emission-line mapping. 

RJS thanks Mike Hudson and Laura Parker for useful 
discussions about this work, and the Anglo- Australian 
Observatory for retrieving the raw AAT spectra of 
ESO325-G004. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. This publica- 
tion makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center /California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. 
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